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ABSTRACT 



It is well-established observationally that the characteristic angular momentum axis on small scales 
around AGN, traced by radio jets and the putative torus, is not well-correlated with the large-scale an- 
gular momentum axis of the host galaxy. In this paper, we show that such misalignments arise naturally 
in high-resolution simulations in which we follow angular momentum transport and inflows from galaxy 
to sub-pc scales near AGN, triggered either during galaxy mergers or by instabilities in isolated disks. 
Sudden misalignments can sometimes be caused by single massive clumps falling into the center slightly 
off-axis, but more generally, they arise even when the gas inflows are smooth and trace only global grav- 
itational instabilities. When several nested, self-gravitating modes are present, the inner ones can precess 
and tumble in the potential of the outer modes. Resonant angular momentum exchange can flip or re-align 
the spin of an inner mode on a short timescale, even without the presence of massive clumps. We therefore 
do not expect that AGN and their host galaxies will be preferentially aligned, nor should the relative align- 
ment be an indicator of the AGN fueling mechanism. We discuss implications of this conclusion for AGN 
feedback and BH spin evolution. The misalignments may mean that even BHs accreting from smooth 
large-scale disks will not be spun up to maximal rotation, and so have more modest radiative efficiencies 
and inefficient jet formation. Even more random orientations/lower spins are possible if there is further, 
un-resolved dumpiness in the gas, and more ordered accretion may occur if the inflow is slower and not 
self-gravitating. 
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1 INTRODUCTION 

Understanding accretion is critical for inferring the origin of the su- 
permassive black hole (BH) population ( Soltan'1 982[|Salucci et al.| 
[1999, Shankaretal. 12004, Hopkins et al. 2006b). Most of the BH 
growth in the Universe is obscured by large columns of gas and 
dust, so knowing the behavior of gas on scales ^ 0.1 — lOOpc is 



also |Natarajan & PringTel|1998) or obscuring AGN "torii" defin- 
ing the plane along which material flows into the inner accretion 
disk (see e.g. |Keel|1980l [Lawrence & Elvisl 1982 1 l Ulvestad & WiT 
"son 1984, Sch mitt et al. |1997", Sim coe et al.|[T 997, Kinney et al. 
2000,,Gallimore et al.|2006t|Zhang et al.|2009f nThe nuclear disk 



a necessary ingredient in a full model of BH evolution ( Antonucci 



T982^ T^93l, "Lawrence 1991] [Risaliti et aLl[T999| [Simpson et ah 
1999, Willott et al. 2000] ). The discovery of tight correlations be 



tween BH mass and host spheroid properties (e.g. mass, velocity 
dispersion, binding energy; "Kormendy & Richstone*1995",'Magor-' 
Irian et a l. 1998, Ferrarese & Merritt 2000, Gebhardt et al. 2000, 
'Hopkins et al. 2007b a. Alter & Richstone 2007, Feol i & Mancini[ 
2009 ) implies that BH growth is coupled to galaxy formation. Mod- 
els widely invoke some form of feedback from AGN to explain the 
origin of the BH-host relations, rapid quenching of star formation in 
bulges, the color-magnitude relation, and the cooling flow problem 
(e.g.|Silk&R ees 1998, King 2003 2005, Di Matteo et al. 20051 
Springel et al., 2005, Hopkins et al. 2008, Hopkins & Elvis,201o| 



is misaligned with the larger- scale disk/galaxy inflows; but the lat- 
ter must ultimately be the origin of the former, so this is not trivially 
expected. 

This misalignment has a number of consequences. It con- 
strains any model of AGN fueling and has important implications 
for AGN obscuration. Not only does it constrain the origin of the 
"torus," but misalignments between the inner and outer disk can 
potentially result in large covering factors of obscuration (even if 
the disks are thin; se e e.g. [ Sanders et al.[[1989[ jNay akshinl[2005 [ 
jFruscione et al.|2005{ [Hopkins et al.|2 01 laV It is critical for under- 
standing BH spin - if gas accreted from large scales in the galaxy 
conserves its axis of angular momentum as it falls onto the BH, then 
almost any high accretion rate event will spin the BH up to near- 
maximum (a ^ 0.998) and align it with the parent disk/inflow (e.g. 



Croton et al.|2006[ and references therein). 

However, despite these important links, the detailed processes 
in BH fueling remain poorly understood. One critical long-standing 
puzzle is the consistent observational finding that there is little or 
no correlation between the angular momentum axis of material ac- 
creting onto the BH, and the axis of the host galaxy. This has been 
seen with a number of different tracers, e.g. radio jets (expected to 
align with the axis of the BH spin or inner accretion disk, but see 
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Volonteri et ari[2005l [Volonteri & Rees'" 2005[ jBerti & Volontenl 



2008j ). However, if the angular momentum can be randomized on 



sufficiently small mass/timescales ("chaotic accretion"), then not 
only will the lack of correlation with the host galaxy appear ( |King[ 
[& Pringlel[2007| ), but the typical spins are held low even in large 



^ We stress that this is not necessarily the same as a lack of correlation 
between obscuration and host galaxy alignment, since significant obscur- 
ing columns can come from large scales in e.g. starbursts or edge-on disks 
jHopkins et al. 2006a Hopkins & Hernquist 2006 Hayw ard et al.|2011a[ 
[Zakamska et al.,2006„Rigby et al.,2006,,Lagos et al.,2011f . 
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accretion events ( [Moderski & Sikora|l 996 1 [King & Pringle|2006| ). 
Spin has important subsequent implications for BH-BH mergers 
and gravitational wave BH recoil (whether or not BHs will be ex- 
pelled from the galaxy or rapidly damp any small recoil motion). 
And it is believed to be critical for the production of radio jets (at 
least in some scenarios; seefBlandford & Znajekf T977l |Begelman| 
let al." 1984 , but also compare Livio et al. 1999 ). Jets and other AGN 
feedback sources are of critical importance for quenching cooling 
in massive galaxies, shaping the galaxy mass function, structuring 
galaxy clusters, and resolving the "cooling flow problem." 

Unfortunately, it is not generally possible to simultaneously 
model inflows from galactic scales and their behavior on the small 
scales near the BH that are relevant for this problem. Analytic mod- 
els ([Kawakatu & Wada|2008l [Nayakshin & King|2007[|Elitzur &| 
|Shlosman|200 6^ are limited by symmetry assumptions as well as 
the fact that these systems are highly non-linear, often chaotic, and 
not in steady- state (with inflow, outflow, star formation, and feed- 
back competing). Simulations of galaxies used to follow inflows 
are typically limited to a resolution of several lOOpc, much larger 
than the scales of interest ( [Cattaneo et al.|[2005l [Hopkins et al.| 
|2005a|b| ). Other simulations which begin on small scales (taking 
some fixed initial conditions for the gas inside of < lOpc), can- 
not relate this to the larger- scale material from which it must have 
originated ( [Schartmann et al.|2009||Wada & Norman|2002| |Wada| 
|et al.|2009] l. Some exciting results have emerged from "zoom-in" 
refinement techniques (see Escala|2007[IColpi et al.|2007[ |Levine| 
let al.|2008[ [Mayer et al.|2007||Dotti et al.|2009| ), but computational 
expense has generally required restrictive assumptions (e.g. turn- 
ing off cooling and star formation on small scales) or limited these 
to single example galaxies at a single instant in time (preventing 
statistical statements). 

Recently, [Hopkins & Quataert| ( |2010a| ) attempted to build on 
these experiments to model the angular momentum transport re- 
quired for massive BH growth, and carried out a series of nu- 
merical simulations of inflow from galactic to BH scales. By "re- 
simulating" the central regions of galaxies in a series of stages, 
gas flows can be modeled over a range of galactic scales from 
^ 100 kpc to < 0.1 pc. In [Hopkins & Quataert[ { |2010 a) we show 
that quasar-level inflows (~ IOMq yr~ ) arise from global pertur- 
bations such as galaxy mergers and/or secular instabilities, which 
(when sufficiently strong) generate a cascade of subsequent in- 
stabilities (of varied morphology), and typically manifest near 
the radius of influence of the BH as a thick (torus-like), lop- 
sided/eccentric gas-i-stellar disk. In [Hopkins & Quataert[ p010b| ) 
we discuss evidence for the rehcs of such disks in nearby galaxies 
(jLauer et al.|1993[ [Bender et al.|2005|[Lauer et al.[T996l ). In [Hop^ 
[kins & Quataert| ( |2011a| ) we discuss the detailed dynamics of these 
instabilities and how they drive large inflow rates, and in [Hopkins] 
|et al. (2011a) discuss their role in the obscuration of AGN. 

In this Letter, we show that the instabilities which drives in- 
flows in these simulations naturally lead to large misalignments of 
the nuclear disk with respect to the disk of the host galaxy. We 
discuss the implications for BH spin, even in "maximally conser- 
vative" scenarios where there is no unresolved sub-grid dumpiness 
in the ISM in our simulations (there almost certainly is such). 

2 THE SIMULATIONS 

The simulations used here are taken from a suite used to study 
the physics of gas inflow from galactic to small scales injHop- 
I kins & Qu ataert ( 2 0T0al [20TTal [20 1 Ob| ) ; [Hopkmsl pOTO] ) ; [Hopkins 
'& Quataert (2011b). The numerical properties of each simulation 
are specifically given in ^Hopkins & Quataert, ( ,201 Oa^ ) (Tables 1- 
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Figure 1. Illustration of twists, warps, and misalignments in some of our 
simulations (each panel is a different simulation). Projected gas density (in- 
tensity) and specific star formation rate (color, increasing from blue through 
yellow) are shown. Times are chosen randomly near the peak of the BH 
accretion. All are projected "face on" to the disk on large scales (angular 
momentum averaged over the entire box). These are ultra high-resolution 
"intermediate-scale" simulations which can resolve twists and misalign- 
ments between the "torus" and larger-scale disk. There is frequently such 
a misalignment or warp between the inflow from larger scale bars-within- 
bars, and the nuclear disk (e.g. middle-left or bottom-left cases), or mis- 
alignment driven by the inflow of large clumps from fragmentation of large- 
scale modes (e.g. top-left or top-center cases). 

3), but we briefly describe them here. In order to probe the very 
large range in spatial and mass scales, we carry out a series of 
"re- simulations." First, we simulate the dynamics on galaxy scales. 
Specifically, we use representative examples of gas-rich galaxy- 
galaxy merger simulations and isolated, moderately bar-unstable 
disk simulations. These are well-resolved down to ^ 100 — 500 pc. 
We use the conditions at these radii (at several times) as the ini- 
tial conditions for intermediate- scale re- simulations of the sub-kpc 
dynamics. In these re- simulations, the smaller volume is simulated 
at higher resolution, allowing us to resolve the subsequent dynam- 
ics down to ^ lOpc scales - these re- simulations approximate the 
nearly instantaneous behavior of the gas on sub-kpc scales in re- 
sponse to the conditions at ^kpc set by galaxy-scale dynamics. We 
then repeat our re- simulation method to follow the dynamics down 
to sub-pc scales where the gas begins to form a standard accretion 
disk. 

Our re- simulations are not intended to provide an exact re- 
alization of the small-scale dynamics of the larger-scale simula- 
tion that motivated the initial conditions of each re- simulation (in 
the manner of particle- splitting or adaptive-mesh refinement tech- 
niques). Rather, our goal is to identify the dominant mechanism(s) 
of angular exchange and transport in galactic nuclei and what pa- 
rameters they depend on. This approach clearly has limitations, es- 
pecially at the outer boundaries of the simulations; however, it also 
has a major advantage. By not requiring the conditions at small 
radii to be uniquely set by a larger- scale "parent" simulation, we 
can run a series of simulations with otherwise identical conditions 
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(on that scale) but systematically vary one parameter (e.g., gas frac- 
tion or ISM model) over a large dynamic range. This allows us to 
identify the physics and galaxy properties that have the biggest ef- 
fect on gas inflow in galactic nuclei. The diversity of behaviors seen 
in the simulations, and desire to marginalize over the uncertain ISM 
physics, makes such a parameter survey critical. 

The simulations were performed with the TreeSPH code 
GADGET-3 (Springel 2005,) ; they include stellar disks, bulges, dark 
matter halos, gas and BHs. For this study, we wish to isolate the 
physics of gas inflow and so do not include explicit models for 
BH feedback (see § |4|. The simulations include gas cooling and 
star formation, with gas forming stars at a rate oc p^^^ moti- 
vated by observations (|Kennicutt|1998| ), and normalized so a MW- 
like galaxy has total M* ^ 1 Mq yr~ . Varying the exact slope or 
normalization of this assumption has no qualitative affect on our 
conclusions. Because we cannot resolve the detailed processes of 
supernovae explosions, stellar winds, and radiative feedback, feed- 
back from stars is modeled with an effective equation of state 
( [Springel & Hernquist|2003| ). In this model, feedback is assumed 
to generate a non-thermal (turbulent, in reality) sound speed that 
depends on the local star formation rate, and thus the gas density; 
the results shown span a wide range in this "effective sound speed" 
without any strong dependence on the exact value (detailed com- 
parisons of the effects on morphology and inflow are shown in |Hop-| 

t> & Quataert 2010a; comparisons of mode growth and torques 
lopkins & Quataert 2011a ). More detailed comparison with the 
explicit stellar feedback models presented in |Hopkins et al.| ( [2bl lb| 
|2012|[2011c| ) will be the subject of future work. 

We "begin" with galaxy- scale simulations that motivate the 
initial conditions chosen for the smaller- scale re- simulation cal- 
culations. These include galaxy-galaxy mergers, and isolated bar- 
(un) stable disks. These simulations have 0.5 x 10^ particles and 



spatial resolution of 50 pc (details in Di Matteo et al. 2005 , Robert- 



son et al.|2006l|Cox et al.|2006[ [Younger et al. 2008 , Hopkins e t al. 
2009| r a subset have ~ 10 particles and 20 pc resolution. From this 



suite we select representative simulations of gas-rich major merg- 
ers of Milky-Way mass galaxies and their isolated bar-unstable ana- 
logues, to provide the basis for our re- simulations. Small variations 
in the orbits or the structural properties of the galaxies will change 
the details of the tidal and bar features on large scales; however, 
we show in ,Hopkins & Quataert| ( t2010a ) that the precise details of 
these large-scale simulations do not instantaneously alter the dy- 
namics on small scales (see Figs. A2 & A3 therein). Rather, the 
local dynamics depends on global parameters such as the gas mass 
channeled into the central region, relative to the pre-existing bulge, 
disk, and black hole mass (set, of course, by the large-scale inflows, 
but once set, robust to variations in the details of that inflow struc- 
ture). 

Following gas down to the BH accretion disk requires much 
higher spatial resolution than is achievable in the galaxy-scale sim- 
ulations. We therefore select snapshots from the galaxy-scale sim- 
ulations at key epochs and isolate the central 0.1 — 1 kpc region 
which contains most of the gas driven in from large scales (typ- 
ical ^ IO^^Mq in gas, over scale-length ^ 0.3 — 0.5kpc). From 
this mass distribution, we then re-populate the gas in the central 
regions at much higher resolution, and simulate the dynamics for 
several local dynamical times. These "intermediate- scale" simula- 
tions involve 10^ particles, with a resolution of a few pc and particle 
masses of^ IO'^Mq. Wehave run ^50 such re- simulations, corre- 
sponding to variations in the global system properties, the model of 
star formation and feedback, and the exact time in the larger- scale 
dynamics at which the re- simulation occurs. ,Hopkins & Quataert, 



P010a| ) present tests of this re- simulation approach and show that 
it is reasonably robust for this problem. This is largely because, for 
gas-rich disky systems, the central ^ 300 pc becomes strongly self- 
gravitating, generating instabilities that dominate the subsequent 
dynamics. 

We repeat our re-simulation process once more, using the cen- 
tral ^ 10 — 30pc of the first re- simulations to initialize a new set of 
"small-scale" simulations. These typically have ^10^ — 10^ parti- 
cles, a spatial resolution of 0.1 pc, and a particle mass ^ 100 Mq. 
We carried out ^ 50 such simulations to test the robustness of our 
conclusions and survey the parameter space of galaxy properties. 
These final re- simulations are evolved for ^10^ years - many dy- 
namical times at 0.1 pc, but very short relative to the dynamical 
times of the larger- scale parent simulations. 

To check that our re- simulation approach has not introduced 
any artificial behavior, we have run a small number of higher res- 
olution "bridging" simulations. These result in slightly worse ulti- 
mate spatial resolution than the net effect of the "re- simulations," 
but they obviate the need for the re- simulation and bridge the 
scales of the above simulation suites. These include 6 simula- 
tions on galaxy scales (3 mergers, 3 isolated disks) with > 10^ 
gas particles and lOpc softening lengths. While not quite as high- 
resolution as our "intermediate- scale" re-simulation runs, these 
provide an important check on the results of the latter and are run 
self-consistently for 4 x 10^ yr. We have followed the same proce- 
dure on small scales: running 5 "intermediate- scale" simulations 
(with a range of gas fraction and bulge-to-disk ratio) with > 10^ 
gas particles and softening of ^ 0.3 pc; thes e extend from scales 
0.3 — 1000 pc and are run for 2 x 10^ yr. In Hopkins & Quataert 



(2010a) and Hopkins & Quataert| pOlla) we explicitly compare 
the results of these simulations with those of our "re- simulations" 
at the dynamic range where they overlap, and find they are very 
similar (see e.g. the discussion and Figs. 9-13 & A4 in "Hopki ns &| 
Quataert 2010a and Fig. 8 in [Hopkins & Quataert|2 01 la), support- 
ing the methodology used for most of our calculations. 

We note that recent studies comparing cosmological sim- 
ulations done with GADGET and the new moving mesh code 
AREPO fSpringer"2010') have called into question the reliability 
of smoothed particle hydrodynamics (SPH) for some problems re- 
lated to galaxy formation in a cosmological context ( Vogelsberger 
et al.|2011[|Sijacki et al.|20TT]|Keres et al.|2011||Bauer & Springel 



2011 



However, we have also performed idealized simulations of 
mergers between individual galaxies and found excellent agree- 
ment between GADGET and AREPO for e.g. gas-inflow rates, star 
formation histories, and the mass in the ensuing starbursts { |Hay-| 
[ward et aL]|2011b |). Simulations of this type circumvent many of 
the issues with SPH by characterizing the gas on small scales with 
an effective equation of state (as in the present study), rather than 
attempting to resolve the various gas phases explicitly. The dis- 
crepancies above are also minimized when the flows of interest 
are supersonic (as opposed to sub- sonic), which is very much the 
case here ( [Kitsionas et al.|2009| [Price & Federrath|2010[ [Bauer &| 
Springel 201 1 ). We have also performed direct resolution studies of 
simulations at each "scale" (with up to 168 times as many particles) 
and find good convergence (see e.g. § Al and Fig. Al in Hopkinsj 
|& Quataert|2010a| and Fig. 4 in [Hopkins & Quataert|2011a| L 

3 RESULTS 

Figure [T] shows the central tens of pc in several of our 
"intermediate- scale" bridging simulations, in which inflows are 
followed from 1 — lOOOpc scales. Of course, resolving those 
larger scales and the resulting inflow means that the resolution on 
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Figure 2. Alignment of the gaseous disks as a function of radii, across our 
sample of simulations (each line is one simulation chosen near the peak in 
inflow). We quantify (mis)alignment as the ratio of j^/ j^, where j is the 
total angular momentum vector of the gas within an annulus around radius 
R, and the z axis is (by definition) the angular momentum axis of the entire 
galaxy. Within < lOpc, there is relatively weak correlation between the 
inner and outer disk angles. 



these scales is not quite as good as our "small-scale" runs, but the 
lOpc scale disk is marginally resolved (in length/mass; the ver- 
tical/internal structure is not resolved below these sizes in these 
runs) There are clear cases where the inflows from sub-kpc scale 
bars map onto the disk at the BH radius of influence, but with a 
very significant misalignment between the two. 

Figures [2] & [4] quantify the degree of misalignment of the nu- 
clear regions in the simulations. Since the observable quantity is 
generally the absolute value of the misalignment, we plot / f 
(i.e. cos^ 0), where j is the specific angular momentum in an annu- 
lus, and the z axis is the axis of the net angular momentum vector 
of the entire galaxy. Figure |2]plots this as a function of radius, at a 
given time in each simulation. Note that there are misalignments at 
all radii from 0.1 pc to lOkpc (although the cases where there 
are significant misalignments on >kpc scales are generally galaxy 
mergers). Figure [4] plots the cumulative distribution of this quan- 
tity at a fixed small radius, summed over all times and the entire 
ensemble of simulations. The misalignments on small scales are 
somewhere between pure random and pure alignment. 

Figure [3] illustrates the time evolution of the inflow axis. We 
plot the evolution of the central angular momentum orientation as a 
function of time: specifically the angle defined between j(t) (total 
angular momentum vector within 5 smoothing lengths of the BH - 
a couple pc) and the (fixed) z-axis (initial j = jz). Variation in cj) 
(azimuthal angle) is much more rapid, but is less significant phys- 
ically (since systems are axisymmetric to lowest order cj) variation 
reflects lopsided/eccentric modes). There is large time- variability. 
The most extreme cases exhibit several "flips" with > 1: /2 (anti- 
alignment of the disk with its original inclination). 

Fully understanding how this affects BH spin would require 



^ See Fig. 10 in |Hopkins & Quataert|2010a| which shows the vertical scale 
heights as a function of radius in these and the others of our simulation 
suite, compared to our SPH smoothing. At ~ lOpc, the disks have scale 
heights from ~ 1 — 3pc, compared to a softening of ~ 0.3 pc, so the internal 
structure cannot be resolved at smaller radii. The true "nuclear scale" re- 
simulations have resolution of ^ 0.1 pc, and so resolve h/R to ~ 1 — 3pc. 




4 6 
Time [Myr] 

Figure 3. Polar angle between the nuclear disk at our smallest resolved 
scale and the initial (uniform) angular momentum axis of the entire sys- 
tem, as a function of time (same simulations as Figure [2] for clarity and to 
show short time-scale variability, we plot only a small fraction of the sim- 
ulated time). The BHs here are accreting at ~ 10% of Eddington; at this 
rate, systems whose inflow angular momentum axis is effectively random 
over ~ 10^ yr timescales (the duration shown here) will be spun down to 
low/modest spin values; this includes most of our simulations. If they ac- 
crete at Eddington, the relevant timescale is ~Myr; this includes only the 
most rapidly variable simulations. 



a number of sub-grid assumptions beyond the model here (see e.g. 
Fanidakis et al. 2011). Even if we assume that the gas retains its 
angular momentum axis below the resolution limit, we need to 
follow the orientation and magnitude of the BH spin as a func- 
tion of time, which evolves as, at first, Lens-Thirring alignment 
forces the inner accretion disk to either align or anti-align (depend- 
ing on whether jdisk • j'bh > or < 0, respectively) inside of some 
warp radius /?warp, and then the torques associated with this even- 
tually re-orient the BH spin in alignment with the disk ( Bardeen &| 
|Petterson]|1975| ). For a given "event," full alignment will occur if 
cos^ > — /j/27bh (where is the original angle between the BH 
and disk, /bh the BH spin angular momentum, and Jd is approx- 
imately the interior disk angular momentum passing through the 
warp region; [Scheuer & Feiler|p^96l [King et al.||2005] ). But this 
also depends on the sub- structure, dynamics, and properties of the 
internal a-disk, well below our best-case resolution (see Kum ar &| 
|Pringle|1985j . However, crudely speaking, for typical a-disk mod- 
els, this translates to a criterion on the mass accreted in a given 
"event" with coherent angular momentum: if the angular momen- 
tum remains coherent over a timescale long enough for the BH to 
accrete some fraction (typically a few percent; |Lodato & Pringle| 
|2006|[Pe"rego et al.[200 9 ) of its mass, then the spin will re-orient to 
align (even if initially retrograde) and most of the accretion will go 
to spinning up the BH. If the inflow angular momentum is incoher- 
ent on this time/mass scale, however, the spin undergoes a random 
walk with decreasing magnitude ( [King & Pring le 2006 , King et aL| 
|2008| ). If the BH is accreting at a fraction A times Eddington, this 
corresponds to a physical timescale of ^ 10^ yr. Consider this 
in Fig. [3] If the accretion is sufficiently rapid {\^ 1), then only the 
most extreme simulated variability will be sufficient to give very 
low spins. However, for the more typical A ^ 0. 1 (coherence time 



10^ yr) in observed systems feollmeier et al.||2006 



Hickox et al. 



[20091 [Hopkins & Hernqui st 2009 , Trump et al.|2009 



and actually 

calculated (via the inflow rate into < 0.1 pc) in these simulations, a 
large fraction of the simulations have sufficient resolved precession 
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Figure 4. Distribution (averaged over time, and over the ensemble of sim- 
ulations) of the alignment/angle (7^ / p = cos^ 9) between the nuclear disk 
and angular momentum axis of the large-scale system (where jz = j, by 
definition). We measure this in three radii: the smallest resolved radii used 
in Figure [3] a fixed physical radius of 0.5 pc, and a radius enclosing ten 
times the gas mass of that enclosed in the smallest resolved radius. The 
three agree reasonably well, suggesting the results are robust at small radii 
(though this clearly reaches the limits of reliable resolution effects). To 
compare, we show the distribution which would be obtained if the orien- 
tations were purely random and uniformly distributed over the sky. Pure 
alignment would be a 5 -function at /j^ = 1 . The distribution is closer to 
random than to pure alignment; still, there is an obvious tendency for more 
alignment than in the pure random case. There may also be a weak excess 
of misalignments near cos^ ~ 0.4 — 0.5, but this is marginally significant. 



in their inflow angular momenta to produce a "random walk" spin 
behavior. 

4 DISCUSSION 

Using high-resolution simulations of gas inflows from galaxy to 
sub-pc scales around AGN, we study the evolution of BH-host 
galaxy alignments. We predict only a weak correlation between 
the nuclear axis and the large-scale disk axis. If anything, this is 
a lower limit to the typical degree of "randomness" in alignment, 
as more clumpy star formation or infall from recycled stellar wind 
material can increase the variation in orientations. Twists and mis- 
alignments, therefore, can explain the random alignment of AGN 
disks relative to their host galaxies. A warped or twisted disk may 
also yield large covering angles towards the BH even when the disk 
itself is thin, although we argue in Ho pkins et al.| ( |20l7a] ) that this 
is not alone sufficient to explain observed obscuration (the "torus" 
must also be geometrically thick). 

These misalignments occur for at least two reasons. First, 
there are cases where large-scale fragmentation occurs in the gas 
(part of a spiral arm or other instability fragments and sinks to 
the center), which can dramatically change the nuclear gas angu- 
lar momentum content (see also |Nayakshin & King|2007[ |Levine| 
|et al.|2010||King et al.|2008| ). And second, even in perfectly smooth 
flows, it is well-known that secondary bars in the presence of dissi- 
pative processes (i.e. gas) will tend to de-couple their angular mo- 
mentum from the primary bar (e.g. [Heller et"aL]|2001| and refer- 
ences therein). Inflow and dissipation lead to runaway strengthen- 
ing of the inner mode, which populates various chaotic orbit fam- 
ilies and exchanges angular momentum with the outer mode, de- 
coupling the inner mode angular momentum and orbit plane from 



that of the outer mode { [Hasan & Norman] 1990[ [Heller & Shlosman] 
T996 , Maciejewski & Sparke|2000| ). The inner mode precesses or 
tumbles in three dimensions relative to the outer mode frame, a 
phenomenon seen in a large number of simulations ( Shlosman & 
Heller|2002[[El-Zant & Sh losman'200'?, 'MacT ejewski & Athanas 



soula[[2008[ [Englmaier & Shlosman 2004 ) and observed double 
(and even triple) bars ( Shaw et al. 1995 , Friedli & Martinet 1993! 



Friedli et aL^1996, ,Erwin & Sparke 1999 , Laine et al. 2002 , Erwin] 
& Sparke|2002i rThese processes are common in our simulations. 



especially in the complicated triaxial potential of realistic merger- 
formed bulges. 

An analogous process also occurs here with the inner lopsided 
disk, at the inner radius (ILR) of the outer bar (itself, in several 
cases, the "inner" of a double bar). [Hopkins & Quataert[ ( [2011a| ) 
show in both these simulations and analytic calculations that an- 
gular momentum exchange in the gas in the central regions (in- 
side the BH radius of influence but outside the viscous accretion 
disk) can be strongly dominated by super- sonic gas shocks sur- 
rounding strong torquing regions in the stellar nuclear disk with 
lopsided/eccentric (m = 1) modes (see also Salow & Static? 
Jacobs & Sellwood|200ll [Sambhus &"Srid ha r 2002. Bacon 



^'2 00T| 
1 et al'[ 



2001[ ). These modes can resonantly exchange angular momentum 



with the pattern at larger radii in the manner of nested bars, leading 
(in plane) to possible reversals and counterrotation of the pattern, 
which in turn reverses the sense of torques on the gas. If the mode 
is strong enough, the exchange in strong shock regions can be large 
enough to change the gas angular momentum by an order-unity 
factor. Generally, as the gas approaches the mode, it experiences a 
sudden, strong resonant torque, shortly followed by or accompany- 
ing a strong shock that dissipates its energy. When the torques are 
sufficiently strong the gas falls in on a nearly radial orbit along the 
pattern; the small "residual" angular momentum can have differ- 
ent signs depending on the instantaneous pattern speed, precession 
rate, and resonance structure of the mode (all of which continuously 
evolve). Moreover, we show in [Hopkins et al.[ ( [2011a| ) that when a 
sufficiently strong m = 1 mode appears, the inner disk becomes 
vulnerable to the "firehose instability," and self-excites large ver- 
tical bending modes. The linear derivation of those modes therein 
suggests that they have both large growth rates and order-unity sat- 
uration amplitudes (i.e. drive order-unity fluctuations in in Fig.js]); 
if there is a significant population of stars in the nucleus on retro- 
grade orbits (from, say, previous accretion episodes), the growth 
rate and saturation amplitude of these modes is greatly enhanced 
( [Sellwood & Merritt,1994>,Davies & Hunter,1997 ). We should note 
that, although the evolution in Fig.[3]is extremely "rapid" relative to 
e.g. secular processes a >kpc radii, it is still indeed secular: at 1 pc 
around a BH of 10^ - 10^ Mq, 1 Myr represents - 100 - 1000 dy- 
namical times, so the relevant resonant effects can collectively op- 
erate over a very large number of orbital periods. All of these pro- 
cesses become more prominent as the nuclear gas is more strongly 
self-gravitating; so they may operate progressively more efficiently 
in higher- accretion rate AGN. They do, however, rely on a com- 
plicated interaction between collisionless and collisional material; 
as such, many will not appear in simulations that do not include 
"live" star formation in the disk (compare e.g. | Escala|2007[ |Colpi[ 
[et al.|2007[ ). At least some observed AGN (with e.g. jets and maser 
mapping of their nuclear regions) appear to have such multiple mis- 
alignments corresponding to structures (nested bars) in their hosts 
(see e.g. [Greenhill & Gwinn|1997} . 

We do not predict perfectly random alignments, as there is still 
some bias towards similar axes obvious in Figure |4] (suggested in 
observations as well, in jBattye & Browne|2009|)Shen et al.|2010t 
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|Lagos_et al.|2011| ). Interestingly, there is also a suggestion in Fig- 
ures [2] & |4| of a preference for misalignments of cos^(^) ^1/2 
{0 = (=b7r/4, zb37r/4)). These relative alignments reflect quasi- 
stable potential surfaces, for example, for an inner gas disk in a 
tumbling prolate quasi- spherical potential; they also form the back- 
bone of "X- shaped" (and some "peanut- shaped") bulges formed by 
bar "buckling" after the vertical motions are pumped by resonances 
in the presence of a nuclear mass concentration or secondary bar, 
like those seen here (|Tohline & Durisen||1982| |Pfenniger|[T984l 
[Pfenniger & Norman] 1990| ). It is not surprising, then, that they form 
the upper envelope for the misalignments seen in the more "quies- 
cent" models here. The much larger mis-alignments seen in Fig. [2] 
and most dramatic "flips" seen in Fig. [3] on the other hand, tend to 
arise from the action of large clumps/fragmentation. 

This can have important implications for the spin evolution of 
BHs. If there are no further twists or clumpy structure beyond what 
is resolved here, the resulting spins will in some cases be maxi- 
mal, but in a large fraction of our simulations would be modest - 
changes in alignment on sufficiently rapid timescale will lead to 
rapid BH growth when the orbits are prograde, but then suddenly 
drops when the orbits are retrograde, producing spins in the range 
|a| ^ 0.1 — 0.9 (see e.g. King & Pringle||2006 i. As these authors 
and others have noted, intermediate spins are interesting because 
they imply modest radiative efficiencies (e ^ 0.05 — 0.2, which 
may be favored by BH luminosity density constraints; 'Wa ng et al.| 
[2009. ) , and reduce by a factor of several the fraction of maxima l re- 
coil "kicks" with Av > 1000km s~^ in major BH-BH mergers |van| 
[Meter et all2010[ [Kesden et aLl|2010| ), although this depends on 
the orientation of the orbital angular momentum which may have 
preferred configurations (e.g. Bogdanovic et al. 2007 ), and is itself 
coupled to the accretion history and feedback efficiency in mergers 
(e.g. [Dotti et al.||20T0[ [Blecha et alT][20TT] ). It has also been sug- 
gested that radio jet power may be modest except at near-maximal 
spins (e.g. |Tchekhovskoy et al. 2009 ). 

If real BHs have very low spins, \a\ < 0.2, some additional 
sub-grid processes are required beyond just what we resolve here. 
Either unresolved sub-grid dumpiness in the ISM that would lead 
to more "chaotic" accretion by increasing the randomness of the 
disk orientations on small mass scales as individual clumps are 
accreted ( [King et al.|2008| ), or further twists/bends/misalignments 
continuing into the a disk (see e.g. [Kondratko et al.|2005| [Green- [ 
[hill et al. [2003] [Kinney et al.[2000| ). There may also be resonant ex- 
changes associated with the pairing process in BH mergers, some 
of which may be important to resolve the "last par sec problem" 
( [Colpi et al.|2007l [Dotti et al.|2009| [Nixon et al.|2011b|a| ). On the 
other hand, if near-maximal spins are the norm, then it suggests that 
inflow is somewhat less random than what we find here; some other 
process, such as slower, more extended accretion from low-density 
diffuse gas which is not gravitationally unstable, may dominate. 

Misalignments can also have dramatic implications for BH 
feedback. An outer disk which is misaligned with the inner disk, 
especially one which has multiple "twists," presents a larger "work- 
ing surface" on which AGN feedback may couple (as opposed to 
an AGN in a single, thin disk). Moreover, many feedback mecha- 
nisms are predicted to have preferential alignments corresponding 
to the spin or nuclear gas disk - radio jets and ionization cones be- 
ing preferentially polar, broad absorption line winds preferentially 
planar. If the inner disk precesses rapidly, these mechanisms might 
appear effectively isotropic to the gas at larger scales in the galaxy. 

The results here are reminiscent of those of IBarnes & Hern^ 
[quist[(T99T] [1996') on somewhat larger scales. They found that the 
angular momentum of gas flowing into the nucleus of a merger 



remnant can lose its memory of the initial direction of the disk an- 
gular momentum]^ On this basis, they argued that kinematically 
decoupled cores in elliptical galaxies may originate in this man- 
ner ([Hernquist & Barnes1[T99T] [Cox et al.[|2006l [Hoffman et al.[ 
[201()| ). In particular, frames from their animated sequences of their 
mergers ( [Barnes & Hernquist[1998| ) often display similarities to the 
images shown in Figure [T] 

The results here represent a first study of inflows in a relatively 
"smooth" medium. In future work, we will extend the models here 
to include the effects of realistic stellar feedback, star formation, 
and ISM structure, as well as more detailed physical models for 
AGN feedback. Stellar feedback should always be present in some 
form, and may (as discussed above) further enhance the "random- 
ness" of the angular momentum on small scales. However, the mi- 
crophysics of star formation and stellar feedback in the vicinity of 
even a quiescent BH, let alone a rapidly accreting QSO, are quite 
uncertain. AGN feedback is potentially important during phases of 
rapid accretion, but this is less clear - it may be that the effective 
duty cycle of strong feedback is such that it is not dynamically im- 
portant for the angular momentum evolution of accreted material 
during the time when most of the mass is actually accreted (as, 
when it becomes strong, it suppresses subsequent accretion). This 
will, of course, depend on the specific feedback mechanisms. We 
also specifically avoid the BH-BH merger stage, choosing to focus 
instead on the more simplified case where there is a single BH in the 
galaxy nucleus. Certainly, ongoing pair merging may introduce ad- 
ditional misalignments (and will have important spin effects, noted 
above), but since misalignments are observed even in quiescent, 
isolated systems, their origin must be more general. 
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